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^ \ We calculate the combinatorial R matrix for all elements of Bi Bi where Bi 

p ^' denotes the ■'-perfect crystal of level Z, and then study the soliton cellular 

D ■ automaton constructed from it. The solitons of length I are identified with 

^ . elements of the A^^^-crystal B^i. The scattering rule for our soliton cellular 

^ automaton is identified with the combinatorial R matrix for A [^''-crystals. 
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I. INTRODUCTION 



A cellular automaton is a dynamical system in which to points in the one- 
dimensional space lattice are assigned discrete values which evolve according to a 
deterministic rule. Soliton cellular automata (SCA) are a kind of cellular automata 
which possess stable configurations analogous to solitons in integrable partial dif- 
\^ ■ ferential equations. Solitons move with constant velocity proportional to the length 
. if there is no collision. After collision their lengths are preserved but the phases are 
shifted. A simple example of SCA is the "box-ball system" ( [l8[ , [l^)) for which 
the scattering rule of solitons was shown to be described by the combinatorial R 
O ' matrix of crystal bases of the quantum group Crystal bases Q) 

are a powerful combinatorial tool in the representation theory of quantum groups, 
which, when they exist, relate the modules of a quantum group to a combinatorial 
structure called a crystal. In 0, a class of SCA was formulated in terms of crystals 
^ I of affine Lie algebras. In this formulation the time evolution operator is given by 
' the row-to- row transfer matrix of integrable vertex models at g = 0, represented 
as a product of the combinatorial R matrix of crystals. In [4], using this approach, 
the scattering rules for the SCA associated with the perfect crystals given in jsf 
were determined. In ^ the phase shift of the U'^{An'') SCA was related to the 
energy function of A^^2^-crystals. In 1^, Yamada gave the scattering rules for the 
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SCA associated with the perfect D4^^-crystals given in 12], and in [igI] those of 
the SCA associated with the Z)!^"*- crystal ^B"'' were also given recently. 

The affine Lie algebra Gg^'' is the infinite dimensional analog of the finite dimen- 
sional exceptional Lie algebra G2 (see For each integer / > 1 it has a perfect 
crystal Bi = which is realized as a set of restricted semistandard tableaux 
with a single row and at most / columns in [20'|, however the explicit action of 
the arrows was not given. In [l5[, a coordinate construction of Bi as a subset of 
(Z/3)^ was given along with the explicit action of the arrows. 

In this paper, using the realizations of Bi given in 2^ and 15 1 we compute the 
combinatorial R matrix Bi^Bi^Bi ®Bi for G^2 \ then use it to find the scattering 
rule of the SCA associated with Bi. We see that solitons of length / are parame- 
terized by the A crystal B^i and that when two solitons collide, the scattering 
rule is described by the combinatorial R matrix B^i^ ® B-^i^^B^i^ ® B^i-^Ji > I2. 
This SCA has two characteristic features. Firstly, the parameterization of length 
/ solitons by B^i is very tricky (Proposition [TOj) . Such a parameterization has not 
been seen before. Secondly, the phase shift of a longer soliton may become negative 
after collision (Theorem H]). This phenomenon is observed only in the D\ -SCA 
(fioj) and ours. 



II. PRELIMINARIES 



Let g = G2 he the affine algebra whose Dynkin diagram is depicted as follows. 

1 2 

Let {ao, ai, ^2} be the set of simple roots, {ho, hi, h2} be the set of simple coroots, 
and {Ao, Ai, A2} be the set of fundamental weights of q. The null root 6 and the 
canonical central element c are given by: 

6 = ao + 2ai + 3a2, c = ho + 2hi + h2. 

Let P = span^jAo, Ai, A2, 5} be the weight lattice. Uq{g) denotes the quantum 
group associated with q. 



A. Crystals 

In this section we give the basic definitions regarding crystals. Our exposition 
follows 0. Let / = {0,1,2}. 

Definition 1. A crystal associated with Uq{Q) is a set B together with maps wt : 
B ^ P, Ci, fi : B ^ B U {0}, and £j, yjj : — t- Z U {—00}, for i E I satisfying the 
following properties: 
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1. (piib) = Eiip) + (/ij, wt{b)) for all i e I, 

2. wt{eib) — wt{b) + ai if e-ib e B, 

3. wtifib) = wt{b) - a, ifhbeB, 

4. £i{eib) = ei{b) - 1, ipi{eib) = (pi{b) + 1 if eib e B, 

5. Eiifib) = £i(6) + 1, ifiifib) = cpiib) - 1 if fib e B, 

6. fib = b' if and only if b = e-ib' for b,b' & B,i G I, 

7. if (fi{b) — —00 for b & B, then e-ib — fib — 0. 

A crystal B can be regarded as a colored, oriented graph by defining 



The tensor product Bi B2 of crystals Bi and B2 is the set Bi x B2 together with 
the following maps: 

1. wt(6i ® ^2) = wt(6i) + wt(62), 

2. ei{bi (g) 62) = max(£i(6i), £1(62) - {hi, wt(6i))), 

3. (pi{bi (g) 62) = max{ipi{b2),(Pi{bi) + (/ij, wt(62))), 



where we write bi <S) &2 for (61, 62) £ -Bi x B2, and understand 61 (g) = (g 62 = 0. 
Bi <S> B2 is a. crystal, as can easily be shown. 

Definition 2. LetBi andB2 be Uq{Q)-crystals. A crystal isomorphism is a bijective 



map * : Bi U {0} B2 U {0} such that 

1. *(0) = 0, 

2. if be Bi and ^(6) G B2, then wt{^{b)) = wt{b),ei{^{b)) = £^(6), = 
(Pi{b) for all i E I, 

3. ifb,b' e e B2 and fib = 6', then fi<^{b) = and = 
Ci'^ib') for all i e /. 



fib = b'. 
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3x3 = (mod 2) 

Xi,Xi,X2 — Xs,X3 — X2 



B. Perfect crystal Bi for G^^^ 

For a positive integer I define the set 

-6;= I 6 = {xi,X2,X3,X3,X2,Xi) G (Z>o/3)^ 

For b — {xi, X2, X3,xs,X2,xi) e Bi we denote 

fh\ J , X3 + X3 X3 + X3 
s{b)=xi+X2-\ \- Xi,t{b) = X2-\ — (1) 

and define 

A = {0, Zi,Zi + Z2, Zi + Z2 + 3Z4, Zi + Z2 + Z3 + 8^4, 2zi + Z2 + Z3 + SZi}. (2) 

Now we define conditions {Ei)-{Eq) and {Fi)-(Fq) as follows. 

( (Fi) Z1+Z2 + Z3 + 3^4 <0,Zi+Z2 + 3Z4 <0,Zi+Z2< 0, Zi < 0, 
(F2) Z1+Z2 + Z3 + 3Z4 <0,Z2 + 3-24 <0,Z2< 0, Zi > 0, 

(F3) 2:1 + 2:3 + 3^:4 < 0, 2:3 + 3^:4 < 0, 2:4 < 0, Z2 > 0, zi + Z2 > 0, , s 

(F4) + ^3 + 3^4 > 0, Z2 + 3^4 > 0, Z4 > 0, ^3 < 0, ^1 + ^3 < 0, ^ 

(F5) 2:1 + Z2 + -23 + 3^4 > 0, Z2 + 3^4 > 0, Z3 >0,Zi< 0, 
(Fq) Zi + Z2 + Z3 + 3Z4 > 0,Zi + Z3 + 3Z4 > 0, ^1 + 2:3 > 0, Zi > 0. 

(Ei), (1 < i < 6) is defined from (Fj) by replacing > (resp. <) with > (rcsp. <). 
Then for b = {xi,X2,X3,X3,X2,xi) G Bi, we define ei{b), fi{b),ei{b),ipi{b),i = 0, 1,2 
as follows. We use the convention: (a)+ = max(a, 0). 
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eo{b) 
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m 
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if Z2 > i-Z3) + , 

if 0:2 < < Z3, 

if {Z2)+ < -Z3, 

if {Z2)+ < -Z3, 
if X2 < < Z3, 

if Z2 > (-2:3)+, 

if Z4 > 0, 
if Z4 < 0, 

if 24 < 0, 
if 2:4 > 0. 



(6) 

(7) 

(8) 
(9) 



X2 - 3^3) + ) + , Vl{b) =Xi + (X3 -X2 + {X2 - ^3) + ) + , 
4., 992(5) = 3X2 + |(X3 - X3) + , 

{2zi + 2:2 + ^3 + 32:4), ^oib) = I - s{b) + max A. 



For b E Bi if b' = Ciib) or b' = fiib) does not belong to namely, if Xj or Xj of 
b' for some j is negative or s{b') exceeds /, we understand b' to be 0. 

Theorem 1 ([3]). For the quantum affine algebra Ug{G^^) the set Bi with the 
maps Cj, fi, Ei, ipi, i = 0,1,2 is a perfect crystal of level I. 



FIG. 1: The crystal graph of Bi. Each ^/ is a 1-arrow, each \ is a 2-arrow, and 

arrows pointing up are 0-arrows. 
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Remark: The ?7g(D|)- crystal is contained in Bi 
Let 



fi = {1,2,21,22,23,3, 0,0, 3,23,22,21,2,1}, 
and define the following partial ordering on B: 

1< 2 < 2i < 22 < ^ ^ ^1 < 22 < 2i < 2 < T. 

Proposition 1 ([io*] 3.9, see also [lo[). The Uq{G2) crystal B{jAi) is given by the 
following set of restricted semistandard Young tableaux whose entries are in B: 

ak < ak+i{k = 1,2, ... J - I) 

t2,{b)+to{b)+t^{b)+t2,{b) < 1 

t2Ab) + t2,{b) + t2,ib) < 1 
tti t2sib)+t2M+t2Ab)<l 
t23(&)+sgn(t3(6))+toW < 1 
t5(6) + sgn(t3(&)) + t23(fe) < 1 J 



where ti = {a^ = i\k = 1,2, . . . , j}L ti = \{ak = i\k = 1,2, . . . , j}L and 
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1 ifn>0, 
z/ n = 0. 



We let denote the empty tableau, and do not put a box around a tableau 
with a single entry. 



Proposition 2 ([2(| prop. 2.1). Forgetting the 0-arrows, Bi = 0^=0 -^(i^i) as a 
Uq{G 2) -crystal. 



The tableau 





^2 







has the coordinate following representation: 



3t2 + 2t2i + ^22 + ^23 + tf. 2t2i + 4)^22 + t23 + 6)^3 + + 3^0 + 3t: 



23 



3 ' 3 

3t23 + ^0 + 3^0 + 6^3 + ^23 + 4t22 + 2t2i k + ^23 + ^22 + 2t2i + 3t2 T 



where tk,tk,to, and ^5 defined above. 



C. Crystal for G^^^ 

In this section we define the [/^((jj^'')- crystal whose crystal graph is given 
below: 








It is the Ug{G2 )-crystal B^'^ corresponding to the vertex 2 of the Dynkin diagram. 

D. Combinatorial R : Bi Bu ^ Bv ® Bi 

In this section, we define the combinatorial R matrix for Bi Bi' and the energy 
function H -.Bi^Bv ^ Z. 

Proposition 3 ([9]). There exist a unique crystal isomorphism TZ : Bi ® Bii ^ 
Bit®Bi and a function H : Bi®Bi' — )■ Z unique up to an additive constant satisfying 
the following property: for any b ^ Bi, b' E Bi', and i E I such that Ciip ® b') ^ 0, 

r H{b®b') + l ift = 0,ifo{b)>eo{b'),Mb')>eoib), 
H{B,{b(^b'))= I H{b®b')-l ift = 0,Mb)<£o{b'),Mb')<eo{b), (10) 
H{b (g) b') otherwise, 

where b' ®b = 7l{b (g) 6'). H is called an energy function on Bi (g) Bu. 

Proposition 4 (Yang-Baxter Equation). The following equation holds on Bi ® 
Bu ® Bu> 

(7^ ® i)(i ® 7^)(7^ ® i) = (i ® 7^)(7^ ® i)(i ® n), (ii) 

where 1 denotes the identity map. 

We define 7^i,i+l to be the map 1®^-^ ® 7^ (g) : (g)Jl^ Bi^ lS)ti ^h- In 

this notation, ffTTl) becomes: 

T^i,i+lT^i+l,i+2T^i,i+l = '^i+l,j+2'^i,i+l'^i+l,i+2 (12) 

for L > 3,z = l,2,...,L-2. 

Definition 3 (Affinization). The affinization of Bi is defined to he the set 

{z''b\n e Z,beBi} 

where the action of Ci, fi is defined as: 

eiiz^'b) = 2"+''«gi(6) 

/,(^"6) = z^-'^^m. 

The combinatorial i?-matrix 7^^^ : AS{Bi) O AS{Bi') Aff(SzO ® Aff(Si) is 
given by: 

7^^^(^"6 ® ^"6') = ® (13) 

where b' ®b = TZ(b ® b'). 
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E. Combinatorial R matrix for A^^^ 
Define Bi to be the set {(xi,X2) G Z>q|xi + X2 = 1} together with the maps: 
eg (xi, X2) = (xi - 1, + 1), e^(a;i, xs) = (xi + 1, Xs - 1), 
/^(xi, X2) = (xi + 1, X2 - 1), /i^(xi, X2) = (xi - 1, X2 + 1), 

(Xi,X2) = Xi+i, V^f (Xi,X2) = X2-i, 2 = 0,1. 

For b E Bi if b' = ef{b) or b' = ff^ib) and some Xj is negative in 6', then we 
understand b' to be 0. This defines a perfect f/^ ( A ^^■')- crystal. 
The crystal isomorphism It : Bi ® Bii ^ Bit ® Bi is given by 

'7^((xi,X2) (g) (1/1,1/2)) = iy'i,y2) ® (a;'i,a;2), 

where 

x'l = Xi + min(?/i, X2) - min(?/2, Xi), (14) 

X2 = X2 + min(?/2, Xi) - min(?/i, X2), (15) 

?/i = ?/i + min(?/2, xi) - min(?/i, X2), (16) 

1/2 = 2/2 + min(i/i, X2) - niin(i/2, Xi). (17) 

To (xi,X2) G Bi we associate the following tableau: 



1 


1 




1 


2 


2 




2 
































3:2 





An energy function H : Bi ® B^ T, \s given by: 

i7((xi, X2) ® (i/i, 2/2)) = - min(?/2, xi). (18) 

The combinatorial /^-matrix : ASi{Bi) ® m{Bv) k?i{Bv) ® Aff(-Bi) for A^^^ 
is given by: 

n^^{z"'b ® Z"6') = ^ (19) 

where b' ®b = 'R{b ® V). 



III. COMBINATORIAL R MATRIX FOR G^^^ 

In this section, we give the combinatorial R matrix for the f/^((j>'2^'')-crystals 
Bi ® Bi and B^ ® Bj. This is done using an algorithm similar to the column 



insertion scheme in 13 . 
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A. Combinatorial R matrix for highest weight elements 



First we find the highest weight elements of Bi ® Bi using the following lemma. 



Lemma 1 ([ll']). An element 61 ® 62 ^ ® ^if^) ^■5 highest weight if and only 
if Cibi = and gf^^'^+^62 = Ofori = l, 2. 



By proposition 2.2, as a Uq{G2) crystal, 



S;®6i = 0i3(jAi)®0i3(A;Ai) 

j=0 fc=0 

I I 
= 0(S(jAi) ® B{0)) © 0(i3(jAi) ® B{Ai)). 

j=0 j=0 



Therefore, the highest weight elements are those of the form (j, 0, 0, 0, 0, 0) b' 
where b' E B{0) or B{Ai) and satisfies ei^''^'^^ {b') = 0. We therefore obtain by 
inspection of the crystal graph in figure 1 that the following is a complete set of 
highest weight elements of Bi® By. 



Highest weight element bi ® 62 


Classical weight 
wt(6i 62) 


(j, 0, 0, 0, 0, 0) ® (1, 0, 0, 0, 0, 0) < J < / 
(j, 0, 0, 0, 0, 0) ® (0, 0, 0, 0, 0, 0) < j < / 
(j, 0, 0, 0, 0, 0) ® (0, 1, 0, 0, 0, 0) 1 < J < / 
(j, 0,0,0,0,0)0(0,1, 1,1,0,0) 1<3<1 
{j, 0, 0, 0, 0, 0) O (0, 0, 1, 1, 0, 0) 1 < J < / 
(j, 0, 0, 0, 0, 0) ® (0, 0, 0, 2, 0, 0) 2<j<l 
(j, 0,0, 0,0, 0)0(0,0,0,0,0,1) 1<J<1 


(j + l)Ai 

(j-l)Ai + 3A2 
(j - l)Ai + 2A2 
jAi 

(j - 2)Ai + 3A2 
(j - l)Ai 



Proposition 5. Let IZ : Bi ® Bi ^ Bi ® Bi he the U^{G^2^) crystal isomorphism 
Then we have, letting bi = (j, 0, 0, 0, 0, 0).- 



7^(&l 



1,0,0,0,0,0 



0,0,0,0,0,0 



1,0,0,0,0,0 
1,0,0,0,0,0 
0,0,0,0,0,0 
1,0,0,0,0,0 



1,0,0,0,0,0 



1,0,0,0,0,0 



1,0,0,0,0,0 
1,0,0,0,0,0 
1,0,0,0,0,0 
1,0,0,0,0,0 



/, 0,0, 0,0,0) z/62 = 
and j 

/, 0,0, 0,0,0) z/62 = 
and j 

J + 1,0, 0,0, 0,1) zfb2-- 

and < 

/- 1,0,0,0,0,0) %fb2-- 
and j 

J, 0,0, 0,0,0) z/62 = 

and < 

/-1,1,0,0,0,0) z/62 = 
and j 

j-1,1,1,1,0,0) zfb2-- 

and 1 < 



J<1, 



j-l,|,i,l,0,0) z/62 

and 1 < 

j-1,0,1,1,0,0) z/62 = (0,0,1,1,0,0 

and 1 < J < 
j-2,1,0,0,0,0) z/62 = (0,0,0,2,0,0 

and 2 < j < I, 
0,0,0,0,0,1) z/ 62 = (0,0,0,0,0,1 

and J = 1, 
J -2,0,0,0,0,0) z/62 = (0,0,0,0,0,1 

and 2 < j < I. 

(3) 



1,0,0,0,0,0 
1,0,0,0,0,0 

/-I, 

1,0,0,0,0,0 
J < ^ - 2, 
0,0,0,0,0,0 

0,0,0,0,0,0 

J<l-l, 

0,1,0,0,0,0 

I, 

0, 1,0,0,0,0 

J<l-h 

0,i,|,l,0,0 



Proof. Recall from the remark following theorem 1 that the -crystal 



(20) 



iS^'' is 



contained in the G2^''-crystal Bi. So when (j, 0, 0, 0, 0, 0) ® 6' is a ^4^'' highest 
weight element, the result follows from proposition 3.7, which leaves us to 
check 

7^((J, 0, 0, 0, 0, 0)®(0, ^, ^, 1, 0, 0)) = (1, 0, 0, 0, 0, 0)®(j-l, ^, ^, 1, 0, 0), 1 < J < /. 

There are two cases to consider. 
Case i: 1 = 1. In this case 

7^((l, 0, . . . , 0) ® (0, ^, ^, 1, 0, 0)) = nfofJ'fiiio, 0, 0, 0, 0, 0) ® (1, 0, 0, 0, 0, 0))), 

= /0/1/2 /l7^((o, 0, 0, 0, 0, 0) ® (1, 0, 0, 0, 0, 0)), 
= /0/1/2 /i(o, 0, 0, 0, 0, 0) ® (1, 0, 0, 0, 0, 0), 

= (1,0,0,0,0,0)0(0,^,^,1,0,0). 



10 



Case ii: I > 1. In this case 

7^((J, 0, . . . , 0) ® (0, ^, ^, 1, 0, 0)) = nif^fJ^Mo, o, o, o, o, o) ® (i, o, o, o, o, o))), 

= Mfiflfimo, 0, 0, 0, 0, 0) ® (1, 0, 0, 0, 0, 0)), 
= ^/i/s 0, 0, 0, 0, 0) ® (1, 0, 0, 0, 0, 1), 

= (l,0,0,0,0,0)®(j-l,i ^,1,0,0). 



□ 



B. Energy function for G2^^-crystal Bi ® Bi 

Theorem 2. The energy function H : Bi ^ Bi ^ Z, I > 1 is given on the highest 
weight vectors by the following: 



H{ij,0,...,0)®b') 



0, if J = I, b' = (1,0,0,0,0,0), 
-1, if j = i_i and b' = (1, 0, 0, 0, 0, 0), 

or j = 1 and b' = (0, 1, 0, 0, 0, 0) or (0, 0, 0, 0, 0, 0), 
—2 otherwise. 

(21) 



Proof. Recall from the remark following theorem 1 that the Z)4^''-crystal B^''' is 
contained in the G2^''-crystal Bi. So when (j, 0, 0, 0, 0, 0) (g) 6' is a D^^ highest 
weight element, we can use the result of corollary 3.8. This leaves us to verify 
that H{{j, 0, 0, 0, 0, 0) ® (0, |, i, 1, 0, 0)) = -2, 1 < j < /. There are two cases to 
consider. 

Case i: 1 = 1. In this case H{{0, 0, 0, 0, 0, 0) ® (1, 0, 0, 0, 0, 0)) = 0, since j = 0. 
Thus, 

H{{0, 0, 0, 0, 0, 0) ® (0, 1, 1, 1, 0, 0)) = /7(/i/|/i((0, 0, 0, 0, 0, 0) ® (1, 0, 0, 0, 0, 0))), 



3' 3' 



i/((0, 0,0, 0,0,0)® (1,0,0,0,0,0)), 

-1. 



Therefore, 



H{il, 0, 0, 0, 0, 0) ® (0, ^, ^, 1, 0, 0)) = H{fo{{0, 0, 0, 0, 0, 0) ® (0, ^, ^, 1, 0, 0))) 
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since /o acts on the first component. 

Case ii: I > 1. Then, H{{0, 0, 0, 0, 0, 0) O (1, 0, 0, 0, 0, 0)) = -2. Therefore 

H{{o, 0, 0, 0, 0, 0) ® (0, ^, ^, 1, 0, 0)) = HChllhm 0, 0, 0, 0, 0) ® (1, 0, 0, 0, 0, 0))), 

= ii'((0, 0,0, 0,0,0)® (1,0,0,0,0,0)), 
= -2. 

By proposition 5 we have 

7^((0, 0, 0, 0, 0, 0) ® (1, 0, 0, 0, O, O)) = (l, O, O, 0, 0, 0) ® (1, 0, 0, 0, 0, 1), 

so by acting on (1, 0, 0, 0, 0, 0) ® (1, 0, 0, 0, 0, 1) by the appropriate sequence of fi 
we see 7^((0, 0, 0, 0, 0, 0) ® (0, |, |, 1, 0, 0)) = (1, 0, 0, 0, 0, 0) ® (0, |, |, 1, 0, 1). Now, 
letting /o act on (0, 0, 0, 0, 0, 0) ® (0, |, \, 1, 0, 0) j times, we get (j, 0, 0, 0, 0, 0) ® 
(0, |, |, 1, 0, 0), where each apphcation of /o acts on the first component. However, 
for the corresponding action on (1, 0, 0, 0, 0, 0) ® (0, |, |, 1, 0, 0), we see that each 
apphcation of /o acts on the second component. Therefore, 

H{{j, 0, 0, 0, 0, 0) ® (0, ^, ^, 1, 0, 0)) = i/((0, 0, 0, 0, 0, 0) ® (0, ^, ^, 1, 0, 0)), 

= -2. 

□ 

C. Combinatorial R matrix for Bi (E> Bi 

Let 6i 62 ® &3 ® ■ ■ ■ ® be a highest weight element of the G2-crystal B{Ki)®^ 
and B{hi 62 • • • h) be the (j2-crystal generated by 61 ® &2 ® &3 ® • ■ ■ ® In analogy 
with the "column insertion" algorithm in [l3[ let ^ : B{1 0) — )■ -B(l),?7 : B{1 12)—)' 
B{1 2 1), and 6 : B{1 1 23) — B{1 23 1) be the unique isomorphisms between the 
given crystals. We remark that the G2-crystal considered in pjj] is not the same 
as our B{Ki)®\ so the analogy is not precise. Let 

: i3(l)^*-i ® B{1 0) ® i3(l)®'-^-i ^ i3(l)®' 
r]i : i3(l)®*-i ® B{1 1 2) ® i5(l)®'-*-2 ^ ^ b{1 2 1)® B{1)®^~'-^ 

Oi : i3(l)®*-i ® i3(l 1 23) ® i5(l)®'-^-2 ^ g(i 23 1) ® i5(i)S5i-i-2 

be the following crystal isomorphisms: 

61^ = 06*® l®^-^-2 
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Lemma 2. Let ai 0:2 ■ ■ ■ be a semistandard G2 tableau, and let /3 G B{1) 
be such that cti ® /5 G 2) (resp. B{1 2^) ). Then the image 771772 • • • iln-i{oin ® 
Oin-i <8) • • • <8) ai (g) ^) ( resp. 6162 ■■ ■ On-\ ® Q;„_i ® ■ ■ ■ ® a^® P)) is well-defined. 



Proof. Since ai 0:2 ■ ■ ■ ctn is semistandard we can apply a suitable sequence of 
e^s so that (8) q;„_i (8) • • • (g) cki (8) ^ becomes l®!®---®!®^' with 1 ^' e B(l 2) 
(resp. 23)). We can then apply a suitable sequence of e^s to get I®!®- • ■®1®2, 
(resp. 1 (g) 1 (g) ■ ■ ■ (g) 1 (g) 23). Now, apply 7711/2 • • • Vn-i (resp. ^1^2 • • • ^n-i) on 
l(g)l®---(g)l®2, (resp. 1 (g) 1 (g) ■ ■ ■ (g> 1 (g) 23). Use the corresponding sequence 
of fiS and commute with 771772 ■■■ 7/„-i (resp. ^^i^^2 • • • ^n-i) to get the image of 

QJn <H) Q;„_i (g) • • • (g) CKi (g) CKi (g) ^. □ 

Now we are ready to prove the main result of this section. 

5 ■ ■ ■ (g) ai E Bi and 1)2 = [i or 



Theorem 3. Let hi 





a2 







aj (g) a,- 



G Bi. Then, the unique Ulj{G 2^) -crystal isomorphism TZ : Bi ® Bi ^ Bi ® Bi is 
given as follows: 



) /3 


«i 


^2 










«2 










012 






1 


«] 


as 


> ■■ ■ 





(g) 61 



7^(6l (g) 62) = < 



Ql+i ® 


gi 


92 




9i 






P<g) gi 






9j 


r 








91 


92 




9i 



pi 



bi (g) 62 



Q2 


(1;-! 







r 


q;2 






«i-i 


9 






91 


92 




9j-i 





ifai®/3eB{l l),j=l, 
%fai®PeB{l l),j=/-l, 
i/ai(g)/3eB(l l),j<l-l, 

ifb2 = 0,J = ^, 
z/62 = 0, j < /, 
z/ai®/3eS(l 2),j =/, 
ly/iere g^+i (g) g/ ■ ■ ■ (g) gi = 
^1^2 ■ ■■Vi-iibi ® ^2), 
ifai®f3e B{1 2),1< J < /, 
where p ® r ® qj ® qj^i ■ ■ ■ ® qi — 

CrV^2 • --Vj-iibi ® ^2), 
z/ai®/5e-B(l 23),j>l, 
u;/iere g^+i ® qj ■ ■ ■ ® qi = 

9i92...9,_iibi®b2), 
ifai®[5e B(l 0), 
a2 ® C(«i ® /9) e B{1 1), and j > 1, 
where p® q — ^'~^{aj),r — ^{a.\ (g) /?), 
ifai®(3E B{1 0), 
^2 (g) ^(tti (g> /3) G B{1 2), and j > 1, 
w/iere p (g) g^.i (g) gj_2 ■ ■ ■ (g> 91 = 

771772 . ..rij-2ij{bi (g) 62), 
i/cci (g) ^ = 1 (g) 1, j > 1, 
m aZZ other cases. 

(22) 
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Proof. By proposition 5, we have: 



7^( V 



( 1® 


1' 








® 


1' 








1 ® 


p'+i 


1 




1 ® 




i 






0® 










1 ® 


1^- 




2 




1 ® 


P- 


1 


2 





1 ® 


P- 


1 


23 




1 ® 


P- 


1 







1 ® 


P- 


'i 


2 




1 ® 


P- 


'1 






I 6i ® 


^2 









if /5 : 
if /3 : 
if /3 : 

if 62 
if 62 

if /3 : 

if /3 : 

if /3 : 

if /3 : 

if /3 : 

if /3 : 



= 0,J < /, 

2,l<j</, 

23, J > 1, 

0, J>1, 
3,J>1 

1, J>1, 



in all other cases, 



where the symbol V denotes that the element h is repeated j times. Thus we 
can easily verify that the statement of this theorem is true for the highest weight 
elements and all that remains is to show that the right-hand side of fl22l) is well- 
defined for all elements of Bi 



ai 


02 







G Bi and bo 



(3 or 



Bi. Let bi 

E Bihe given. If G B{1 0), then ^{ai®P) is well-defined by the definition 

of ^. If 02 ® ^{ai ® /5) e B{1 2), and j > 1, then 771772 . . . r]j_2^j{aj ® aj-i ® ■ ■ ■ 
ai = ?7i?72 . . . ?7j_2(Q!j <8)aj_i ® ■ ■ ■ (g)Q;2 ® /?)) is well-defined by lemma 2. 

If e 2) (resp. B{1 23)) then r/1772 • • • 77^-1(0^ (g)aj„i O ■ ■ ■ Oai (resp. 

6162 ■ ■ ■ 6j-i{aj ® Oj-i ® ■ ■ ■ ® ai ® /3) is well-defined by lemma 2. The remaining 
cases are clearly well-defined since they do not use any of the maps ^, 77, or 6'. □ 



Example: Let 61 = |2|23|2i| , 62 = 0, / = 4. Then ai ® /3 = 2 ® G 2) and 
j < I. We compute: 

^i^mi^i ® 23 ® 2 ® 0) = ^i\{2i ® 23 ® 3 ® 1) 

= er^(23®3®3® 1) 
= 2i(8)2(g)3(8)3(8)l 



Therefore 7^(6l O 62) = 2i O 1 3 3 2 
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D. Combinatorial R matrix for Bu (8) Bi 



Recall the t/g(G'2^'')-crystal which is isomorphic as a [/5(G2)-crystal to B{A2). 
Then, we have the following: 

Proposition 6. The complete set of highest weight vectors and their weights is 
given below: 



Highest weight vector in 


Classical 
weight 


Highest weight vector in 

Bi®B^ 


Classical 
weight 




Ai + A2 




Ai + A2 




2A2 


1®3' 


2A2 


1' (8) 


A2 


(8) 1' 


A2 


1'® 


A2 


10 0' 


A2 



Proof. As a C/g(G'2)-crystal B^®Bi^ B{K2) ® (i3(Ai) e S(0)). By lemma 1, the 
highest weight elements in ® Bi are those of the form b®b' where b is highest 
weight in i3(A2) and b' e B{Ki),B{Q) satisfies ef'^^b' = 0. Similarly, the highest 
weight elements in Bi ® B^^ are those of the form 6 6' where b is highest weight 
in B{Ai),B{0) and b' e ^(As) satisfies 



e-'" ■ 'b' 
lib' = 



if 6= 1, 
if & = 0, 



for 2 = 1, 2. The proposition follows by inspection of the crystal graphs of B\ and 
B^. □ 

Proposition 7. Let : B\^® B\ ^ B\ ® B^^ be the unique crystal isomorphism. 
Then TZ is given on the highest weight elements by: 



7^(l' b') 



Proof. By property (2) of the definition of crystal isomorphism, we have the fol- 
lowing: iib®b' e B^i^Bi and ^(6® 6') e then wt (^(6® 6')) = wt(6®6'). 
Therefore, since 1' (8) 1 and 1 (8 1' are the only elements of B\^ (8) Bi and Bi (8 B\^ of 
weight Ai + A2, we must have 1Z{1' (8) 1) = 1 (8 1' and, similarly TZ{1' (8) 2i) = 1 (8 3'. 
For TZ{1' <8) 0), we compute 




g^/ogigo(l'® 0) 
n{elfoeMl' ® 0)) 
elfoheo(n{l'®0)) 
n{l'®0)) 



l'®2i 
^(l'®2i) 

1®3' 

/o/ieo/|(l<8)3') = 0<8)l'. 
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Finally 7^(1' (8) 6) = 1 ® 0', because it is the last remaining highest weight element. 

□ 



IV. SOLITON CELLULAR AUTOMATON 

We define = {&i ® &2 ® ■ ■ ■ ® &l e Bf^\bn = 1, for n sufficiently large} to be 
the set of states of the G^^ sohton cellular automaton. We depict the operation 
n{b ®b')^b'®b by: 



b' 



b' 

Fix Z > and let ui e Bi be the element {I, 0, 0, 0, 0, 0). For p = 6i (8) 62 <8) • • • <8) &l e 
Vl we define the time evolution operator Ti{p): 

Tiip) ®ui = Ull+iR-l-il- • -7^23^12 (m; ®p)- (23) 

The transition of phase (8) &2 <8) • • • <8) ^l) = &i <E) &2 <8) • • • ® &l is depicted as 



bi 



u 



(1) 



= ui. 



bi 



We define the state energy to be the sum: 

L-l 

Ei{p)^-J2H{u^'^^bi+i). 



i=0 



(24) 



We can also use the combinatorial it!-matrix TZ : B\^ (S> Bi ^ Bi <S> to define an 
operator Tt, similar to Ti, by 



m?) ® b{p) = 7^LL+l7^L-l L ■ ■ ■ 7^237^l2(l' ® p). 



(25) 



In this case, b{p) is dependent on the state p so we indicate this dependence in the 
definition. 
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A. G2 -solitons and their scattering rules 

Experience from other sohton cellular automata has shown that states p G Vl 
satisfying Ei {p) = 1 correspond to the so-called "one-soliton states" . 

Proposition 8. In the G'^^ soliton, Ei{p) = 1 if and only if p ^ 1'^^ and is one 
of the following: 

• 1®' ® 3^™ ® 2i ® 2®" ® 1®'^-^ 

• 1®' ® 3^™ ® 22 ® 2®" ® 

/or some 1, 171,72, r G Z>o such that l + m + n + r = L. 

Proof. Let p = hi®h2®- ■ -^bi be such that Ei{p) = 1. For b®b' G i3i(g)Si it is the 
case that n{b®b') = b®b' . Thus, Ei{p) = -H{l^bi) -J2LiH{bi^bi+i) = 1. We 
compute i/(l(g)l) = 0,if(l®2) = i7(l(g)2i) = i7(l(g)22) = i/(l(g)3) = H{1®0) = 
— 1 and if(l(g)6) = —2 in every other case. Therefore, the first 7^ 1 must be either 
2, 2i, 22, 3, or 0. However, H{0 (g) 6) > for all 6 G Si, which rules out that case. 
Similarly, we have H{3 (g) 1) = H{3 (g) 2) = if (3 ® 2i) = if (3 (g) 22) = H{3 (g) 3) = 0, 
or H{3 (g) 6) > otherwise. Therefore, only 1, 2, 2i, 22, or 3 can follow 3. In turn, 
we see that only 2 or 1 can follow 2i or 22, and only 2 or 1 can follow 2. Finally, 
a trailing 1 can only be followed by another 1, or Ei{p) would be greater than 1. 
Therefore, we see that if Ei{p) = 1 then p 7^ 1®-^ and falls into one of the three 
cases given in the statement of the theorem. The converse is easy to prove. □ 

Proposition 9. Letp G Vl- If (iiij)) 7^ then e2Ti{p) = Ti{e2{p)) and Ei{e2{p)) = 
Ei{p), otherwise e2T;(p) = 0. The same relations hold for f2- 

Proof. It is easy to see that £2(^0 = V^2(^i) = 0, so by the tensor product rule 
f2{ui ®p) = ui 0J2P and f2{Ti{p) ®ui) = f2Ti{p) ® ui. Therefore, acting on both 
sides of fl2^ by /2 and commuting /2 with TZ gives: 

f2Tl{p) ®Ui = TZl L+lV-L-l L ■ ■ ■ ^23^12 (Wi ® hp) 
^{Ti{f2p)®Ui, ifflpy^O, 

\ 0, otherwise. 

Therefore, we see that f2Ti{p) = Ti{f2p) if /2P 7^ and f2Ti{p) = otherwise. 
Exactly the same argument works for the 62 case (using the tensor product rule 
for 62 instead of the one for /2). 
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Let p = 6i(g)62®- ■ -^bi, Po = z%i(^z%2<^- ■■<^z^hL G Aff(i3i)®'^, andp' = f2P- 
Then 



'T-.Aff 'T-.Aff 



■ntMi{z%i^p',) 



Mi- 



On the other hand, 



)Aff 



= 7^f^+l7^t!, , ■ ■ ■ nint^{z\, ® hp,) 



Comparing the powers of z in fl2B]) and flTTI) we see that 



L-l 



L-1 



i=0 



i=0 



(26) 



(27) 



therefore, 

£;K/2(p)) = Eiip). 

A similar argument shows that Ei{e2p) = Ei{p). □ 

Recall the A crystal Bi = {{xi,X2) G Z>o\xi + X2 = I}, where {xi,X2) can be 
associated with the tableau: 



1 


1 




1 


2 


2 




2 








XI 








X2 





Proposition 10. Let b G B31 and define a sequence (61, 62, • • • , 
i/ie following triples to the corresponding elements of Bi : 



hi) G B\ by sending 



1 


1 


1 


^ 2 


1 


1 


2 


^ 2i 


1 


2 


2 


H->22 


2 


2 


2 



3. 



T/ien t/ie map ii : i^a/ — )■ Bf given by ii{b) = bi ® ® 
relations 

ii{e^b) = e2ii{b), iiCfib) = f2ii{b). 



bi satisfies the 
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Sketch of proof. In Bi 

2 4 2i A 22 4 3. 
InBs 



1 


1 


1 


1 


1 


1 


2 


1 

— > 


1 


2 


2 


1 

— > 


2 


2 


2 



Now use the definition of /2 on tlie / fold tensor product to obtain the result. □ 

Remark: The above map is a bijection from B31 to {p e Vl\Ei{p) — 1} by 
proposition 9. 

A state of the following form is called and m-soliton state: 

...[h] [l2]...---...[lm]... 

where h > I2 > ■ ■ ■ > Im, ■ ■ ■ ■ denotes a local configuration of the form 

3®m (g, 2®'-™, 3®™ ® 2i 2®'-"^-\ or 3®"^ (g) 22 (g) 2®'-"^-\ m > 0, and the [k\ are 
separated by sufficiently many I's. 

Proposition 11. Letp he a one-soliton state of length I. Then 

1. Ek{p) = min(A;,/), 

2. Tkip) is obtained by rightward shift by min(A;, /) lattice steps. 

Proof. By applying 62 sufficiently many times, p becomes 2®' (8) 1 • • • (8) 1. By 
theorems 2 and 3 we obtain: 





® 2) 


= 1 ® 


p- 


~1 


2^-1+1 


if i > 0, 


(28) 


n{¥ 


® 2) 


= 2® 


2' 








(29) 


7^( r 2'-' 


® 1) 


= 2® 




2i-»-i 


if i < Z, 


(30) 


n{¥ 


(8)1) 


= 1 (8) 


1' 








(31) 



and 



H{ V 


2 


(— j ^ 


8) 2) = -1 if i > 0, 


(32) 




2' c 


3 2) = 0, 


(33) 


H{ r 


2 




8 1) = 0. 


(34) 



where, as before, the symbol V means that b is repeated j times. \i k <l then 

Tk{p) ®ui = nLL+i--- 7^237^12(M/ ® 2®' ® 1^-') 

= l®'^ ® 2®' ® 1^-*^-' ® Mfc, 

and 

L-l 
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Otherwise, 



and 



Tkip) ®Ui=nLL+l--- 7^23^12(M^ ® 2®' ® 1 



I L-L\ 



L-1 



Ei{p) = -Y,H{u^'^®hi)=l. 



1=0 

The result follows because /2 commutes with and preserves Ei (proposition 
9). □ 

We now consider the two-soliton case 

p=...[h]...[k]... 

where h > I2. We can use proposition 10 to associate a two-solitori state T^*(p) := Pt 
at time t with the element z~^^hi ® z~''^b2 G AS{B3i-^) ® AS{B3i^), where ki : = 
— min(r, Zj)t + 7j, where 7^ is the number of 'I's to the left of [li]. If r > I2, then 
we expect to see the longer soliton catch up with and collide with the shorter one, 
and, after sufficiently many time steps separate out into another two soliton state. 
At the end of this section, we will prove that is the case. In the following lemmas, 
we identify elements in AS{B3i^) (g) Aff(i33;2) with their corresponding two-soliton 
states. 



Lemma 3 (Analogous to lemma 4.15 in [19|). Suppose we have a one-soliton state 
p = ...[/].. . corresponding to z~^{i^ 31 — i) G Aff^Bsi), < i < 31. Then T^^{p) is 
another one-soliton state, corresponding to z~''{3l,0), if i = 31, and z^'^^^{i + 
1,3/ — 1) otherwise. We also have b{p) = V if i = 31 and b{p) = 2' otherwise. 

Proof. We ignore the initial 'I's. If i = 3/ then p = 2^'- ® i^^-'-. We compute T^(p) 
as follows: 



1' 



1'- 



1'- 



1'- 



1' 



v. 



Therefore T^{p) = 2®' ® 1®^"' = ^-^=(3/, 0) and b{p) = 1' as required. 

If z = 3r, 1 < r < Z then p = 3®'-*^ ® 2®'' ® 1®^"'. We compute T^(p) as follows: 



1' 



3'- 



3' 



■3' 



3' 



3'-^ 3' •■■3' 
22 2 



3'- 



2'- 



2' ■■•2' 



2' 
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Therefore T^(p) = 1 O 3^^''"'^ ® 22 ® 2®'" O 1^-'-^ = z'''~\3r + 1, 3/ - 3r - 1) and 
b{p) = 2' as required. The cases where r = 0,1 are similar. 

Ifz = 3r + l,l<r<Z then p = 3®'-^'-i O 22 ® 2^'' ® l^^'^^ We compute T^(p) 
as follows: 



1 



1 



■3 — 3' 
3 



■2' 



Therefore T^(p) = 1 ® 3®'" 



2i : 

. 2i ® 2®^ ® 1 



L-l-l 



-k-1 



(3r + 2, 3/ - 3r - 2) and 



b{p) = 2' as required. The cases where r = 0, / are similar. 

If i = 3r + 2, 1 < r < / then p = S®^-""^ ® 2i ® 2®^ ® 1^-^-^. We compute T^(p) 
as follows: 



1 



r 



3' 



3'---3' 



3' 



Therefore Tt|(p) 



3 3 

gCg)/— r— 1 



3'— 3' 
2 

r>, o(g)r+l A 



3'- ■■3' 



3' 



,-fc-i 



2'- 



2' ■■■2' 



2' 



(3r + 3, 3/ — 3r — 3) and 



b(p) = 2' as required. The cases where r = 0, / are similar. 



□ 



Now consider the two soliton state ^^^^^(3/1,0) (g) z^'^'^{i,3l2 — i). Thanks to 
lemma 3, the action of has no effect on the first soliton, and therefore we have 
the following corollary. 



Corollary 1 (Analogous to lemma 4.15 in [19[). Suppose we have a two-soliton 



state p = . . .[h 



[^2] . . . corresponding to z ^(3/i,0) 



Affi^sh) ® Ajf{Bsi2),0 < i < 3/2- Then Tt,(j9) is another two-soliton state, cor- 
responding to ^-^=1(3/1,0) O ^-'=2(3/2,0), z/z = 3^2, and z-''^{3li,0) z-^^-\i + 
1,3/2 — 1) otherwise. We also have h{p) = V if i = 3li and h{p) = 2! otherwise. 

Lemma 4 (Analogous to lemma 4.17 in [19]). Assume that li > I2. For p = 

^-^1(3/1, 0) ® z-''^{i, 3/2 -i) e AjJiBsi,) ® AffiBsi^), 0<i<3l2 we have 

1. mn^ffip)) =n^S{T^{p)), and 

2. b{p) = b{n^ff{p)), 

where we use the shifted energy function H = 2I2 + H in IZ^^ . 



Proof. Case 1: i = 0. Then p 



y—kl I 



(3/1, 0) ® z~''^{3l2, 0). Therefore, 



T^in'^^ip)) 



z 



(3/2,0) ® 2 
(3/2,0) ® 2- 



(3/1,0), 

(3/1,0) (by corollary 1), 
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and b{n^^{p)) = v. 
On the other hand, 

= ^"'U3/i,0)®^-'^(3/2,0), 

and b{p) = V. Therefore T^{n^^{p)) = 7^^*^(T^(p)), which is (1), and b{p) = 
b{n^^{p)), which is (2). 

Case 2: <i< Sk- Then 

n^^{p) = z-''^+^'^-\3l2, 0) ® z-'''~^'''+\3h 
mn^^ip)) = z-^^+^^''~\3l2, 0) ® z-^^-^^''+'-\3li - z + 1, z - 1) (by corollary 1), 

and b{1Z^^{p)) = 2'. On the other hand, 

T^(p) = 2-^=1(3/1,0) z-'''-\3l2 - i + l,i - 1) (by corollary 1), 
n^^{T^{p)) = z-''^+^^^-\3l2, 0) ® z-^^-^^^+'-^{3h - i + 1, « - 1), 

and 6(p) = 2'. Therefore T^(7?.^^(p)) = 7t^*^(T^(p)), which is (1), and b{p) = 
b{n^^{p)), which is (2). □ 



Lemma 5 (Analogous to lemma 4.18 in 19|). Let p G Vl, I > 0, L ^ 0. Then 

1. T^{Ti{p))=Ti{T^{p)), 

2. b{Ti{p))®ui = n{ui®b{p)). 
Proof. By definition 

T^(TKp)) ® b{Ti{p)) = nL,L+l ■ ■ -^23^12(1' ® Ti{p)), 

and 

Tlip) 0Ui= TiL,L+l ■ ■ ■ "7^23^12 (m« ® P)- 
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Now, combining the definitions, we compute: 

T^{Ti{p))^b{'I\{p)^^ui 

= 1^1,1+1 ■ ■ ■ 7^23^127^L+l,L+2 • • • 7^34'^23(1' ^ Ui ® p) 
= 7^L,L+l7^L+l,L+2 ■ ■ ■ 7^23^34^12^23(l' ® Ui (g) p) 

7^347^23'7^347^23^^23'7^127^23^12^(1' ® ® p) 

(by the Yang-Baxter equation) 

= 'R'L+l,L+2'R-L,L+l'R-L+l,L+2 " " " 7?.23^347^127?.23 ® l' ® 

(since TZ ""^(1' (g) ui) = ui® 1') 

= ^L+l,L+2^L,L+l ■ ■ ■ ^23'^12'7^L+l,L+2 " " " ^34^23 («« ® 1' ® p) 
= 7^L+1,L+27^L,L+1 ■ ■ ■ 7^23^12 (m« ® T[,(p) (g) 6(p)) 
= ^L+i,i+2(T/_(T^(p)) ®Ui® b{p)) 

= Ti{T^{p))®Tl{ui®b{p)). 

Comparing both sides we see that T\^{Ti{p)) = Ti(T\^{p)) and b(T^^{p)) ®ui = lZ{ui® 
b{p)). □ 

Theorem 4 (Analogous to theorem 4.9 in [l9[). Let p = . . . [li] . . . [I2] ■ ■ ■ E Vl be 
a two-soliton state with li > I2, corresponding to z^^bi ®z^'^b2 G Aff{Bi^) ® AffiBi^) 
where ki,k2 < 0. Then after sufficiently many applications ofTr,r > I2 the new 
state is given by 

n^S{z^^bi ® Z^^b2) = Z^'^b2 ® z^'^bi 

with phase shift 

k'^-k2 = ki-k[ = 2I2 + H{bi ® 62). 

Proof. Let 

p = z^^{x,,X2) ® z^^{y,,y2) G Aff(i33iJ ® Aff(^3iJ. 

By proposition 9, commutes with ei,ff. Thus, it is enough to check the 
scattering rule for the highest weight elements 

We will show the statement is true by induction on ^2- 

Suppose y2 = 0. Then yi = 8/2- In this case, the cellular automaton time evolution 
is reduced to that of the A^^^ cellular automaton by fl28|) - fl3T|) : 



7^( V 2'-* ® 2) 


= 1 ® 




-1 


2/-i+l 


if i > 




= 2® 


2' 








7^( p 2'-' ® 1) 


= 2® 




2i-i-i 


if i < Z, 




= 1 ® 


1' 
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Therefore, the scattering rule is the same as for the A\'' case, namely: 



T;(p)=7t^^(p), 

for sufficiently large t. 

Assume the statement is true for < ?/2- By corollary 1, T\^{p) = z'^^(3Zi,0) (g) 
z^^~^{yi + 1,1/2 — 1), so the inductive assumption holds. Therefore, T^*(Tt|(p)) is a 
2-soliton state, for sufficiently large t, and 

T*(T^(p))=7^nT^(p)). 

Therefore, by lemmas 4 (1) and 5 (1), we have 

T^(T*(p))=T^(7^nP))• 

Also, by lemma 5 (2), we have h{lZ^^{jp)) = h{jp), and by corollary 1 h{j>) = 2' since 
p is assumed to be a 2-soliton state and 1/2 > 0. By lemma 4 (2), b(Tr{p)) ® Mr = 

TZiUr® h{p)) = n{Ur ® 2'). But e2{Ur ® 2') = Ur 1', TZ{Ur ® 1') = 1' ® Mr, 

and /2(1' ® Mr) = 2' ® Wr- Therefore b{Tr{p)) = 2' and, by repeated application 
of lemma 4.2 (2), we have b{T^{p)) = 2'. Therefore, bCR^^ip)) = b{T^{p)) = 2', 
and T^{Tl{p)) = T^(TZ^^{p)). Therefore, the crystal isomorphisms in can be 
inverted to get 

T^ip) = n'^^p) 

is a 2-soliton state. 

□ 

Corollary 2. Let p = . . . [li] . . . [I2] [Im] ■ ■ ■ be an m-soliton state with 

li > I2 > ■ ■ ■ > Im- Then the scattering of solitons is factorized into two-body 
scattering. 

B. Examples of scattering of Gg^^ -solitons 

In this section, we give some examples of two and three soliton scattering gen- 
erated by computer. Here we take / = 10 in Ti, and L = 55 in Vl- 
Example 1: 

t = 0: 3222211112i2111111111111111111111111111111111111111111111 
t = l: 111132222112i21111111111111111111111111111111111111111111 
t = 2: 11111111322212i221111111111111111111111111111111111111111 
t = 3: 11111111111322ll2i222111111111111111111111111111111111111 
t = 4 : 1111111111111322llll2i22211111111111111111111111111111111 
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The initial state corresponds to 2;°(7, 5) ® z ^(5,1). After interacting, the two- 
soHton state that emerges is 5) ® z~'^{\\, 1) = Td^^^z^i^, 5) ® ^"^(5, 1)), in 

agreement with theorem 4. 
Example 2: 

t = 0: 332ill2i21122lllllllllllllllllllllllllllllllllllllllllllll 
t = l: 111332il2i2122llllllllllllllllllllllllllllllllllllllllllll 
t = 2: IIIIII322I3223IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII 
t = 3: 11111111322ll232i21111111111111111111111111111111111111111 
i = 4 : 1111111111322l2il32i21111111111111111111111111111111111111 
i = 5 : 1111111111113222illl32i21111111111111111111111111111111111 
t = 6 : 111111111111112232illll32i21111111111111111111111111111111 
t = 7 : 11111111111111122l32illlll32i21111111111111111111111111111 
t = 8: 111111111111111122ll32illllll32i21111111111111111111111111 

The initial state corresponds to -2°(2, 7) ® 2;~^(5, 1) ® z~^{l, 2) and the final, sepa- 
rated, state corresponds to ^"^(1, 2)0^-^(2, 4)®z-3(8)(5, 4) = n^^n^in^^{z°{2, 7)® 
2;~^(5, 1) (H) z~^{l, 2)), in agreement with corollary 2. 

We now give an example which shows that it is necessary to have h > h > h 
in corollary 2. 
Example 3: 

t = 0: 2222111133112211111111111111111111111111111111111111111 
t=l: 1111222211331122111111111111111111111111111111111111111 
i = 2 : 1111111122223311221111111111111111111111111111111111111 
i = 3 : 1111111111112233112211111111111111111111111111111111111 

t = 4 : 1111111111111111331122111111111111111111111111111111111 
t = 5 : IIIIIIIIIIIIIIIIIIITII2222IIIIIIIIIIIIIIIIIIIIIIIIIIIII 
t = 6: IIIIIIIIIIIIIIIIIIIIITIIII2222IIIIIIIIIIIIIIIIIIIIIIIII 
t^7 : IIIIIIIIIIIIIIIIIIIIIIITIIIIII2222IIIIIIIIIIIIIIIIIIIII 

Remark: Although separation does not occur in the previous example, if we were 
to apply the result of corollary 2 to this example we would expect to get the 
following: 

n^in^^iz'^iu, 0) ® z~\o, 6) ® z-^\Q, 0)) 

= ^zf (^"'°(6, 0) ^'(6, 6) ® z-i'(6, 0)) 
= z-^^{6, 0) ^-^(0, 6) ^"^(12, 0). 

The reason for the observed lack of separation is perhaps due to the fact that the 
leftmost two solitons in the initial state are shifted oppositely to each other, thus 
creating an overlapping state. 
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